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Abstract 

A  theoretical  derivation  of  the  characteristic  of  a  spherical  probe  in  a  flowing 
Maxwellian  plasma  is  given  together  with  various  limiting  forms  of  the  equations. 
The  limitations  imposed  on  the  validity  of  the  theory  by  the  characteristics  of  the 
sheath  are  discussed.  The  theory  is  then  extended  to  the  case  of  an  n-grid  sensor. 
Finally,  the  application  of  the  theory  to  the  determination  of  plasma  parameters 
from  the  probe  characteristic  is  discussed. 
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An  Analysis  of  the  Behavior  of  a  Multi-Grid  Spherical 
Sensor  in  a  Drifting  Maxwellian  Plasma 


I.  INTRODUCTION 

The  use  of  a  probe  to  determine  the  parameters  of  a  plasma  has  long  been  one 
of  the  standard  methods  of  experimental  measurement.  The  foundations  of  probe 
theory  were  laid  down  in  a  classic  series  of  papers  by  Langmuir. 

In  this  paper,  the  details  of  the  calculation  are  given  for  the  current  flow  to  a 
sphere  in  a  plasma  containing  charged  particles,  where  these  particles  have  a 
Maxwellian  velocity  distribution  with  a  superimposed  drift.  The  current  is  obtained 
as  a  function  of  the  potential  difference  between  the  sphere  and  the  plasma  and, 
following  Langmuir,  the  retarding  and  accelerating  cases  are  treated  separately. 
The  discussion  in  Section  7  imposes  an  important  limitation  on  the  conditions  under 
which  the  expressions  can  be  considered  valid.  The  general  theory  of  an  n-grid 
sensor  is  developed  in  Section  10,  and  the  particular  case  of  a  two -element  sensor 
is  discussed  in  Section  11.  This  last  section,  using  various  special  cases  ol  the 
expressions  developed  earlier,  discusses  the  problem  of  obtaining  the  parameters 
describing  the  plasma;  namely —density,  temperature,  mass,  and  plasma  potential, 
from  an  observed  current -voltage  curve.  A  summary  of  the  equations  with  the 
numerical  constants  most  useful  in  practical  applications  is  given  in  the  appendix. 


(Received  for  publication  11  December  1068) 
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2.  THE  INTEGRAL  FOR  THE  CURRENT 

The  current  crossing  an  element  of  area  dS  is  given  by  J  •  dS  where  J  is  the 
current  density.  If  J  is  due  to  the  motion  of  charged  particles  then  we  car^write 
J  s  Nev ,  where  N  is  the  charged  particle  density,  e  the  charge  on  each  of  the 
particles,  and  v  their  velocity.  If  the  charged  particles  have  a  velocity  distribution, 
then  the  number  of  particles  with  velocities  in  the  range  v,  v  +  dv  may  be  written 
N  f(v)  dv,  and  the  current  density  due  to  this  group  of  particles  is  Nef(v)  dv.  The 
current  flowing  across  the  element  of  area  dS  is  Nef(v)  dv(v  •  dS).  The  total  current 
flowing  across  a  surface  S  due  to  all  of  the  charged  particles  can  be  written 

I  = 


where  the  surface  integral  is  taken  over  the  surface  under  consideration,  and  the  velocity 
integral,  over  the  range  of  velocities  possessed  by  particles  which  cross  that  surface. 

We  consider  a  spherical  collector  of  radius  R  situated  in  a  gas  in  which  the  velocity 
distribution  of  the  charged  particles  is  Maxwellian  and  there  is,  in  addition,  a  super¬ 
imposed  steady  drift  motion  relative  to  the  sphere.  To  calculate  the  current  to  the 
sphere,  we  surround  it  by  a  concentric  spherical  surface  of  radius  r  and  calculate  the 
current  flowing  between  the  two  spheres.  The  problem  now  consists  of  developing  the 
appropriate  distribution  function  for  a  Maxwell  distribution  with  a  superimposed  drift, 
calculating  the  appropriate  limits  for  the  velocity  integral  in  terms  of  the  potential 
difference  between  the  two  surfaces,  and  finally  evaluating  the  integral. 


If 


Ne  f(v)  dv(v  •  dS) 


3.  The  distribution  function 

The  Maxwell  distribution  function  may  be  written 
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where  v'  ,  v'  ,  v* ,  are  the  components  of  the  velocity  parallel  to  a  rectangular  set 
x  y  z 

of  axes,  x'y'z',  at  rest  in  the  gas;  and  a  *  m/2kT  where  m  is  the  mass  of  the  particle, 
T  is  the  temperature  of  the  gas,  and  k  is  Boltzmann's  constant.  If  we  now  super¬ 
impose  a  steady  drift  Vq  in  the  z  direction,  then  the  velocities  vx,  v  ,  vz,  in  a  set 
of  axes  xyz  at  rest  in  the  sphere,  become  vx  =  v'x,  vy  s  v*  ,  and  vz  s  v'z+  vQ. 

To  calculate  the  c’ -Trent  crossing  a  spherical  surface  of  radius  r  we  require 

o 

v  •  dS  which  can  be  written  vrr  sin  Od0d<|>,  where  vr  is  the  radial  component  of  the 
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velocity,  and  0  and  $  are  the  colatitude  and  longitude  angles  in  conventional  spherical 
coordinate  notation.  We  must,  therefore,  express  the  distribution  function  in  terms 
of  the  radial  velocity.  The  velocity  of  the  particles  can  be  represented  in  terms 
of  vr>  vQ,  v  which  are  the  components  in  a  local  cartesian  system  representing 
the  rectangular  components  of  velocity  in  the  directions  of  increasing  r,  0 ,  and  $  at 
the  particular  point  concerned.  Since  vx,  v  ,  vz  and  vr,  vQ  ,  v^  are  two  sets  of 
rectangular  components  relative  to  axes  which  are  at  rest  in  the  sphere,  then 


2  .  2  .  2  _  2  2 
vx  +  vy  +  vz  "  vr  +  Vc 


'0 


and 


'  2  +  v'  2  +  v'  2  =  V  2  +  v  2  +  (v  +  vn) 
x  y  z  x  y  z  0 


*vx+vy2  +  vz2  +  v02+2v0vz 


2  2  2  2 

*  vr  +  vQ  +  +  v0  +  2vq  (vr  cos  0  -  vQ  sin  0 ) 


In  the  transformed  system,  the  volume  element  dvxdv^dvz  becomes  dvrdvQ  dv^ 
since  this  is  still  a  rectangular  system.  To  evaluate  the  integral,  a  further  transform¬ 
ation  is  required  into  a  spherical  coordinate  system  in  velocity  space.  The  colatitude 
e  is  measured  from  the  vr  direction  and  the  longitude  A  is  the  plane  defined  by  vQ 
and  v^  from  the  vQ  direction. 

Carrying  through  the  appropriate  substitutions  gives 


2  2  2  2  2 

v'x  +  v'  +  v'z  =  v  +  vQ  +  2vv0  (cos  c  cos  0  -  sine  cos  A  sin0) 

where  v  is  the  total  velocity  of  the  particle. 

o 

The  volume  element  now  becomes  v  sin  e  dv  de  dA  and  the  integral  for  the 
current  crossing  the  sphere  is 


2  *  jt 
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I  -  f  f  f  f  exP  f-ajv2  +  Vq  +  2vvq(cos  e  cos0  - 

*^>s  0  *0  =0  v  e  A  1 

-  sin  A  cos  A  sin0)jj(vcos  e)(v2sine)  dAdedv  (r2  sin0)  d0  d<J> . 
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The  convention  followed  so  far  is  that  the  positive  direction  is  outward  from  the 
center  of  the  sphere,  so  that  the  current  I  given  by  this  expression  i3  the  current 
flowing  outward.  Since  we  are  interested  in  the  current  flowing  to  a  spherical 
electrode  in  the  plasma  we  need  the  current  flowing  inwards  and  it  is  convenient  to 
change  from  the  convention  at  this  point.  The  velocity  v  is  always  positive  but  we 
now  measure  its  colatitude  from  the  inward  directed  radius  vector,  rather  than  the 
outward.  Hence,  if  we  write  e1  =  nr  -  £,  then  c'  varies  over  the  range  0  to  nr/ 2  for 
inward  flowing  currents.  The  expression  for  I  now  becomes 


2  .  2 
v  +vQ 


-  2vvq  (cos  e'  cos<t>  + 


+  sin  e1  cos  X  sin  0 


v*3 (-cos  e')  sin  e1  dX  (-de1)  dv 


We  shall,  henceforth,  drop  the  prime  from  £,  it  now  being  understood  that  e  is 
measured  from  the  inward  directed  radius  vector. 

The  integrand  is  not  a  function  of  <J>,  which  we  would  expect  from  the  axial 
symmetry  of  the  problem  about  the  drift  velocity.  Hence,  we  can  write 


3 

2 

I  =  27rr2Ne^j 

( cos  e  cos  0  +  s in  X  cos  X  s  in  0  )  J 

Note  that  the  symmetry  of  the  problem  also  requires  that  the  same  current  be 
obtained  regardless  of  the  sign  of  Vg.  (This  is  not  immediately  obvious  from  the 
above  expression. ) 

We  now  need  the  limits  for  the  velocity  integral. 


j  cos  e  sin  e  sin  0  dX  de  dv  dO 


exp(-av2)  J *  J*  j "  Jv3  exp  |-a  |v2  -  2vvg 


0=0  v 


I.  MMITS  OK  IM'KtrK  \TIO\ 

We  assume  that  the  velocity  distribution  of  the  particles  at  radius  r  is  given  by 
the  distribution  function  derived  above.  We  also  assume  a  potential  difference  V 
between  this  surface  and  the  spherical  collector  radius  R.  We  now  calculate  the 
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possible  range  of  velocities  a  particle  can  have  at  radius  r  if  it  is  to  be  able  to  reach 
the  collator  radius  R. 

The  calculation  divides  into  two  cases:  the  retarding  case  when  the  sign  of  the 
charge  on  the  particle  and  the  sign  of  the  potential  on  the  sphere  are  the  same,  and  the 
accelerating  case  when  they  are  opposite.  These  two  cases  were  distinguished  by 
Langmuir. 

Considering  first  the  retarding  case,  then  a  particle  with  velocity  v  travelling  at 
angle  c  to  the  (inward) radius  vector  at  radius  r  will  arrive  at  the  collector  radius 
R  with  velocity  v^  travelling  at  an  angle  to  the  radius  vector  where  these  quantities 
are  related  by  the  laws  of  conservation  of  angular  momentum  and  energy.  If  the 
potential  difference  between  the  collector  and  the  surface  radius  r  is  V,  then  from 
energy  considerations 


|mv2  =  jmv2  +  |eV| 

1  2  I  I 

We  define  a  velocity  v^  by  the  relation  tj-  mvp  =  |eV|  and  this  equation  then 
becomes 

2  2  2 
V  =  V,2  +  VR 

From  angular  momentum  considerations 
rv  sin  e  =  Rvj  sin 

The  particle  will  be  collected  if  it  arrives  at  the  collector  travelling 
tangentially  (sin  ^  ^  1),  that  is,  if 

rv  sine  <  .. 

— H - V1  ' 

which  on  substituting  in  the  energy  relation  gives 


2  2  2 >  r2v2  sin2c 

v  R  ‘  1  ^2 
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i 


\ 


}  I 


Hence,  for  a  given  value  of  v,  the  particle  is  collected  if 

•ta2«a  7  • 

with  the  limiting  value,  sin  being  given  by  the  equality  sign.  The  limits  for  the 
€  integration  are,  therefore,  0  and  €p. 

Note  that  the  upper  limit  of  c  is  a  function  of  v,  one  of  the  other  variables  in  the 

o  . 

integration.  We  obtain  the  limits  for  v  by  noting  that  the  smallest  value  of  sin  c  is 
2  2 

zero  and,  hence,  v  -v^>Oorv>  vR.  The  limits  on  v,  therefore,  become  Vp 
and oo  .  Finally,  for  the  limits  on  X  we  note  that  the  aijuments  given  above  are 
valid  for  all  values  of  X,  and,  hence,  the  limits  for  the  X  integration  are  0  and  2n. 
The  expression  for  the  current  now  becomes 


exp  |-a  jv2  -  2wq(cos  e  cos  0  +  sin  e  cosXsin0)J  j  cos  esincsin0  dXd^dvdO  . 

For  the  accelerating  case,  e  and  V now  have  opposite  signs  and  the  energy  relation 
can  be  written 

o 

jmv2  *  ^mv2  "  le^  ! 

1  2  l  l 

Defining  a  velocity  v'^  by  the  relation  2'mvA  =  leV|»  we  have 


Following  the  same  steps  as  before, 


2  x  ,  2>  r2v2sin2€ 

v  '  V  »  - n - 
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The  integral  must  now  be  evaluated  in  two  parts.  The  above  relation  is  satisfied 
for  all  values  of  c  from  0  to  *72  if 


2  ,2 

v  +  v'A  >  — FT- 


a  r2  -  R' 


,2  „  2  „2 


n  n 

Defining  at  =  R*/(r*  -  R*)  and  VA  =  a  v'A  ,  the  limits  for  the  first  part  of 

A  A  1  o 


integral  become  0  to  *72  for  e  and  0  to  v^  for  v.  (Note  ^-niv^  =  a  |eV|.)  The 
current  I  may  now  be  written 


/  /A  JT  ^ 


I1=27rr2Ne(|)  exp(-a v$)  f  f  j  f  v3 

^0=0  J0  Je=oA=0 

exp  |-a|v2  '  2wq(cos  e  cos  0  +  sin  e  cos  X  sin  0  )jjcos  e  sin  c  sinO  dXdedv< 
If  the  particle  has  a  velocity  greater  than  v^,  it  will  be  collected  if 


•  2  <R2  /,+V a 

Sm  C—  — it  I  1  +  — rr 


^oTTr 


Again,  the  limiting  value  is  given  by  the  equality  sign.  Hence,  the  second 
part  of  the  integral  becomes 
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.*  «  €a(v)  y.  2?r 


I2  =  2»r2Ne^|)  exp(-av2)  f  f  f  f  ^ 

*0=0  vAJ  0  X=0 

exp  |-a|v^  "  2Wg  (cos  €  cos  0  +  sin  e  cos  X  sin  0  )|j  cos  e  sin  e  sin  0  dXde  dvd0  . 


5.  EVALUATION  OF  THE  INTEGRAL 

The  three  integrals  formulated  in  the  last  section  are  identical  as  far  as  the  0 
and  X  integrations  are  concerned,  so  we  evaluate  these  first.  The  integration 
required  is 


it  2 » 

I  I  exp  2a  wQ ( cos  e  cos  0  +  3in  e  cos  X  sin0)  sin0  dXd0 
0=0  a=0  *  J 


Figure  4.  Limits  of  Integration 
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We  can  regard  this  as  an  integration  over  the  unit  sphere  in  0X  space.  The 
direction  cosines  of  a  point  on  the  unit  sphere  are  given  by  /  =  sin0  cos  X, 
m  =  sin0  sinX,  n  =  cos0.  Writing  A  =  2a  wQ,  we  have  for  the  integral 


j^£xp  A  (n  cos  c  +  /sin  e)  dS, 


where  S  represents  the  unit  sphere  and  dS  is  the  element  of  area  sin©  dXd0  in 
spherical  polar  coordinates.  Making  a  cyclical  notation  of  the  direction  cosine  gives 


j^Cxp  A(n  cos  c  +  /sin  e)  dS  =  fj exp  A  (/cos  e  +  m  sin  e)  dS  = 
s  s 

ir 

=  J  J  exp  |  A  (sin  0  cosX  cos  e  +  sin©  sinX  sinc)j  sin©  d0  dX  = 
0=0  x=o 

fn 

-  J  J  exP  JA  sin0  cos(X  -  e)  Jsin0  d0  dX 


0=o  x=o 


Substituting  X'  s  X  -  c  gives 


ff 

0=0  X'  =  -< 


exp  A  sin©  cos  X'  }sin0  d0  dX’ 


and,  since  the  integrand  is  periodic  in  X1  and  we  are  integrating  over  a  complete 
period,  any  complete  period  can  be  used.  Therefore, 


fir  -2 ir-€  ^*2* 

J  J  expjAsinO  cosX'Jsin©  d0  dX’  =  J  J  expj  A  sin©  cosX'  | 
0=0  X'  =  -€  0=0  X'=0 

sin©  d0  dX1  =  ft  expjAsinO  cosxj  sin©  d9  dX  =  ff exp(A l)  dS 

a  =  n 


sin©  d0  dX'  = 


0=0  x=o 
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=  jj  exp  (An)  dS  =  /  /  exP  (A  cos  0  )  sin0  d0  dX  = 

s  0  =0  X=0 

=  exp  (A  cos  0  )  Jo  =  -  Jexp  (-2avvQ)  -  exp  (2avvQ) 


Substituting  the  above  result  into  the  expression  for  current  gives 
3 

7  /* 

I  =  27rr2  Ne  exp  (-av  2)  J  Jexp  (-av2)  J  exp  (-2aw0)  -  exp  (2awQ)J 


V  € 


cos  e  sin  e  de  dv 


and  so 


I  =  27rr2  Ne  l J  J  v2  jexp  [  -a  (v  -  vQ)2  J  -  exp  J  -a(v  +  V())2j| 

0  v  c  ' 

cos  e  sin  f  de  dv  . 

We  now  continue  this  integration  for  the  retarding  case.  For  the  e  integral 
we  have 


f 


1  .  2 

cos  e  sin  e  de  =  j  sin 


=  1  R2  /  VR  \ 

7  7  V  '7/ 


Substituting  this  result  gives 


I  =  v  R  Ne 


oo 

VI  \  J  (v2  -  VR  >  [expj-a(v  -  v/ {  -  exp  J-a(v  +  y/jj  dv  . 


These  integrals  may  be  evaluated  in  terms  of  the  error  function  which  we  define  as 
as 

x 


erf  (x)  exp  (-t2)  dt  . 
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!  ! 
j  5 


;  ? 
i  l 


Using  this  definition,  we  have 


/  y2  expj  -a<v  -  v0>2  |dv  +vo)l  |erf  a?  (VR  ‘  v0)|  +  2T(vR  +  v0) 

exp  j-a(vR  -  v0)2{  ( 

/  y2  [e*p  j-a<v  -  v/j  -  exp  |-a(v  +  vQ)2|  ]dv  =  j  (35  +  vq] 


[' 


1  1 

T  T 

erf  a  (vR  +  vQ)  -  erf  a  (vR 


-  v0>]  +  exP  [-a(vR  -  vo)2J- 


vr  ■  vn  r  2 

-TiT2-  Cxp-a(vR  +  v0) 


/  VR  exP  |-a<v  -  v0)2| dv  =  1  VTv R  t1  -  erf  a?(vR  -  v0>] 


J  VR  exp  j-a(v  -  v0)2J  -  exp  j-a(v  +  vQ)2  Jj  d\  * 

VR 

=  \  VI  VR2  [erf  *  (VR  +  v0>  -  erf  a?  (vR  -  V0>]  • 

On  substituting  in  these  results  and  simplifying,  we  get 


I=xR2Ne^f  (4+v0-vR2)  jerfa2 


<vR  +  v„)  -  erf  a  <vR  -  v„)  }  + 


+  -^Ta-2-  exP  |-a(vR  '  v0>2 1  -  exP  I  ‘a(vR  +  V*) 


V„  -  Vr 
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It  is  convenient  to  define  two  new  parameters 


2  2  m  2 

y  =  avo  ■  2ET  v0 


Finally,  we  have 


I  .  *R2  Ne^  I  ^  +  r2  -  x2  )  [erf  (x  +  V)  -  erf  (X  -  y)  J  + 

+  ^<x  +  y)  exp  { *(x  -  y)2  j  (  X  -  y)  exp  j-(x  +  y>2  j  • 

Turning  now  to  the  accelerating  form  of  the  expression,  the  «  integral  for  the 
first  part  is  simply 


cos  e  sin  e  de  =  -j 


giving  for  1^ 


^  .  jrr2  Ne  i  j*  v2  [exp  |-a<v  -  v0>2}  -  exp  {-a(v  ♦  v/}]  dv  . 


This  can  again  be  evaluated  in  terms  of  the  error  function 

0>>  -  ‘  VI fe  +  Vo)  [6r£  ^  <VA  -  V  +  6rf  Jv0 


VA 

f  v2  exp  j  -a(v  -  vQ, 

*'o 


1  1 


+  J®  exp  <-av2>  -  ^  a.—  exP  { 'a<vA  '  v0) 


and 


14 


f  v2  [exp{-a(v  -  v/  j  -exp  j -a (v  +  v„>2  jj  dv  =  \  VFfei  *  v2) 


1  1  1  1  v 

erf  a7(vA  -  vQ)  -  erf  a7  (vA  +  vQ)  +  2  erf  a7  v0  +  IT  exp  (*avo *  ' 


VA  +  V0  exp!.a(vA-vn)2UVA'V° 


“2T 


A  -V  |+  25  exp  |'a(vA  + V  j 


!-a(vA+vn)2!  . 


For  the  contribution  of  the  first  part  of  the  integral  to  the  current,  this  gives 


xl  =  "2  Ne  iff  ^  [  7  4I(£  +  vo2  ){ert  a*  <v 


1 

A  -v0>  -erfa7<vA  +  v0)  + 


1 

7 


2V  VA+V0 


+  2  erf  a  vQj  -»  —  exp  (-aV())  -  — ^ -  exp  j-a(vA  -  vQ)  J  + 

+  -A25-°  exp{  'a(vA  +  V2}]  ■ 


The  c  integral  for  the  second  part  is 


/ 


cos  €  sin  c  de 


.  1  .2  .1  a  (.  VA  \ 

-7sin  CA  ~  2  Z-TT  \1  +  ^2) 


giving  for  Ig 

I2  *  *r2  Ne  i-  jj-fy  ^ j  [exp  J-a(v  -  V())2j  -  exp  j-a(v  +  v0)2{ 

VA 

This  is  identical  to  the  expression  involved  in  the  retarding  form  of  the 
expression,  except  for  the  values  of  the  constants.  Hence,  we  have  immediately, 


dv. 
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i 

? 

r 

I 

! 


It  is  also  worth  noting  that  the  coefficient  outside  the  whole  expression  contains  the 
radius  of  the  collector  R  r.  the  retarding  case,  and  the  radius  of  the  outer  spherical 
surface  r  in  the  accelerating  case. 


6.  LIMITING  FORM  OF  THE  EXPRESSIONS  FOR  NO-ORIFT 

The  expression  for  no -drift  c»n  be  obtained  from  the  equations  derived  here  by 
taking  the  limit  vQ  —  0,  that  is,  y  —  0.  This  forms  an  important  cneck  on  the 
expressions,  as  the  equations  for  no -drift  are  relatively  easily  derived  from  first 
principles. 

Considering  first  the  retarding  expression,  we  require 


lim  ij'y  y2  *  x2j|erf(x  +  y)  -  erf  (x  -7  )  j  +  ^  (x  +  y)  expj  -(x  -  y)2J  - 

(X  -  Y)  exp  |-(x  +  y)2j]  -  £“,[*  (  ‘  +  r2  -  x2)-erf(x+^-erf(x  '  ?>  + 
+  ^x^P{-(x  -T)V  exp  {  -  (x  +  y)2}  +lsxp  j.(x.y)2j4  lexpj.(x+y)2 

=V^(?  "  *2)jf  GXP  ^"X2)  +  ^  X  4X  ^  ^"X2)  =  2  eXP  ^_X2) 

Hence  on  substituting,  we  obtain 


I  =  2ttR2  Ne^L  exp  ( -x2)  . 


The  accelerating  expression  can  be  dealt  with  in  a  similar  manner  and  involves 
the  same  limiting  expressions 

iff  n .  « -*  r  ™ i-v!!  -TT-rffU*'2--2) 

|erf  (x  +  y)  -  erf  (x  -  y)J  +  ^  (x  +  y)  exp  |-(x  -  y)2|  - 
(X  -  r)  exp  j-(x  ♦  V)2  jj]= 2£r+  exP  <->2)  - 


1 

r 
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1  \J*I  1  „2\  erf(x  +  y)  -  erf(x  -  y)  .  1  exp[  -U  -y)?)  -  exp(  -(x  -  y)2]  . 

(TTT|  )  y  T  y 

+  ^ exp  |-(x  -  y)2  ]  +  ^  e*P  ['<x  +  V>2]  |]  *  2  -  cTTr[f  (?  '  x2)j?  exp  ('x2)  + 

+  2x2  exp  (-x2)  +  exp  (-x2)  J  =  2  |  1  -  --^-y  exp  (-x2) 

Hence,  for  the  current  we  have 

I  =  2*r2  ^  J1  -mexP(-x2)j  • 

Note  that  for  x  =  0  (that  is,  V  =  0)  both  expressions  reduce  to  the  same  equation, 
namely, 

I  =  2ttR2  Ne  ^  . 

^rar 

7.  THE  SHEATH 

Referring  back  to  the  calculations  we  have  made,  we  see  that  in  setting  up  our 
original  integral  we  considered  particles  crossing  a  spherical  surface  of  radius  r. 

In  defining  the  limits  of  integration  we  determined  which  of  these  particles  is  capable 
of  reaching  a  concentric  spherical  surface  of  radius  R  where  R  <  r  when  there  is  a 
potential  difference  V  between  the  two  surfaces.  We  also  considered  only  inward¬ 
travelling  particles  in  setting  up  the  limits  of  integration.  This  means  that  we 
preclude  the  possibility  of  the  potential  on  the  inner  sphere  influencing  an  outward - 
travelling  particle  at  radius  r  in  such  a  way  as  to  reverse  its  outward  motion  so 
that  it  can  eventually  reach  the  sphere  radius  R.  This  is  equivalent  to  supposing 
that  the  electric  field  due  to  the  potential  on  the  inner  sphere  does  not  extend  beyond 
a  distance  r.  This,  then,  is  very  close  to  Langmuir's  concept  of  the  sheath  radius. 

In  deriving  the  above  expressions,  we  have  assumed  that  we  can  equate  the 
potential  to  zero  over  the  sphere  radius  r,  that  is,  that  no  particle  outside  this 
spherical  surface  can  experience  any  force  from  the  influence  of  the  field.  We  have 
assumed  that  the  velocity  distribution  of  the  inward -travelling  particles  at  radius 
r  follows  Maxwell's  law.  (The  velocity  distribution  of  the  outward -travelling 
particles  is  immaterial,  as  they  do  not  enter  into  the  calculation. )  In  deriving  the 
expressions  for  the  limits,  we  used  the  energy  and  momentum  of  the  particle  at 
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the  two  surfaces  only  and  applied  the  conservation  laws.  We,  therefore,  implicitly 
assumed  that  there  are  no  particle  collisions  between  the  two  surfaces.  Also  in 
applying  these  results,  it  is  assumed  that  the  potential  distribution  between  the  two 
surfaces  is  such  that  the  radial  velocity  of  the  particle  nowhere  becomes  zero,  that 
is,  that  there  is  no  potential  barrier  which  the  particle  cannot  surmount.  This  is 
worth  considering  in  a  little  more  detail  as  it  has  important  applications  in  considering 
a  multi-grid  sensor. 

Hie  argument  follows  closely  that  used  to  obtain  the  limits  of  the  integral. 

Suppose  we  have  a  particle  at  the  sheath  edge  at  radius  r.  travelling  with  velocity 
Vj  at  an  angle  to  the  radius  vector,  and  suppose  that  at  some  arbitrary  value  of  r, 

R  <  r  <  r.,  the  velocity  and  angle  are  v  and  e.  For  the  particle  to  reach  the  collector, 
not  only  must  sin€c  (the  value  at  the  collector)  be  less  than  unity,  but  sin  e  must  be 
less  than  unity  for  all  values  of  r  within  the  sheath.  This  means,  from  the  arguments 
given  earlier,  that 


This  condition  will  apply  in  both  retarding  and  accelerating  cases  if  e  and  V 
are  considered  to  contain  the  sign  as  well  as  the  magnitude  of  charge  and  voltage. 
But  the  relation  between  the  limiting  values  of  ei  and  v if  the  particle  is  to  be 
collected  at  all,  satisfy  the  equation 


where  V£  is  the  potential  of  the  collector.  Hence,  eliminating  v. 
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Hie  most  stringent  case  to  be  satisfied  occurs  for  the  accelerating  condition 
when  cj,  can  assume  the  value  x/2.  Hence,  the  condition  is  satisfied  for  all  cases 
if 


which  on  rearranging  gives 


Since  Vc  and  V  may  be  either  positive  or  negative,  this  result  expressed  as  a 
condition  for  V  becomes 


Note  that,  since  we  defined  the  potential  in  the  plasma  outside  the  sheath  to  be 
zerQand  also  R  <  r  <  r.,  the  limiting  condition  implies  that  0  <  |v|  <  |vj  . 

8.  POKER  SERIES  EXPANSION  OF  THE  ACCELERATING  EQUATION 

A  useful  expression  can  be  obtained  by  expanding  the  expression  for  the  acceler¬ 
ating  case  as  a  Taylor  series  in  powers  of  V,  valid  near  the  origin  V  =  0.  The  first 
three  terms  in  the  series  are 


where  the  subscript  0  means  the  value  of  the  quantity  concerned  at  V  =  0.  The 
expression  derived  earlier  gives  I  as  a  function  of  x  rather  than  of  V.  A  compli¬ 
cation  arises  now,  since  we  defined  x  using  the  magnitude  of  eV  — 
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It  is  convenient  to  calculate  the  derivative  (dl/dV)  as  (dI/dx)(dx/dV),  but  in 
differentiating  the  equation  which  relates  x  to  V,  the  sign  of  dx/dV  is  ambiguous. 
This  is  readily  resolved  by  noting  that  the  sign  of  the  derivative  is  the  same  as  the 
sign  of  e,  that  is,  it  is  positive  for  positively  charged  particles,  and  negative  for 
negatively  charged  particles.  The  following  expressions  are  now  quite  straight¬ 
forward  to  obtain 

IQ  =  jtR2  Ne^=g-  I  (1  f  2y2)  erf  y  +  y  exp  (-y2)j  . 

=  tjtR2  Ne  i  2  erf  y  .  Positive  particles 

0  "P*  y  Negative  particles 

Hence,  the  expansion  for  I  becomes 


I  =  *r2  Neja  i[f('  +  2r2> 


erf  y  -  y  exp  (-y  )t 


2|erf *!($)-?“  Tiexp('y2)l@ 


or  with  some  rearrangement  of  the  constants, 


I 


7tR  Ne  v 


|(‘*57) 


7y  exp  (-y2)+  Jerf  V 


12  11  ,  2.|/eV\2  1 

■?V^y!exp('r  MM  •  •]  ■ 

On  referring  back  to  the  complete  expression  for  the  accelerating  case  we  note 
that  it  contains  the  "sheath  radius"  r  rather  than  the  collector  radius  R.  We  also 
recall  that  r  is  the  radius  of  the  spherical  surface  on  which  we  can  assume  V  =  0 
and  this  will  presumably  be  a  function  of  V,  the  voltage  on  the  collector.  Hence, 
we  do  not  have  I  explicitly  expressed  as  a  function  of  V  for  the  accelerating  case 
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until  we  have  r  as  a  function  of  V.  However,  in  the  series  expansion,  the  collector 
radius  R  appears  outside  the  expression  and  r,  which  occurs  in  a,  does  not  appear 
until  the  term  in  V  .  Therefore,  the  linear  approximation  (the  first  two  terms)  is 
independent  of  the  behavior  of  the  sheath.  Unfortunately,  we  cannot  estimate  how 
large  the  voltage  must  be  before  the  linear  approximation  breaks  down,  since  the 
next  higher  order  term  does  depend  upon  the  behavior  of  the  sheath. 

It  is  readily  verified  from  the  retarding  expression  that  the  values  of  Iq  and 
(dl/dV)Q  are  identical  to  those  obtained  from  the  accelerating  expression.  Thus, 
the  curve  of  I  versus  V  is  continuous  and  so  is  its  derivative  at  V  =  0,  that  is,  as 

we  go  from  retarding  to  accelerating.  There  is,  however,  a  discontinuity  in  the 
2  2 

second  derivative  (d  I/dV  )g. 

•».  LIMITING  FORM  OF  THE  EXPRESSIONS  FOR  LARGE  DRIFT  VELOCITIES 

The  behavior  of  the  expressions  for  large  drift  velocities  can  be  obtained  from 
the  power  series  expansion.  This  expansion  was  derived  above  for  the  accelerating 
case,  and  since  Vg  =  yMT,  this  may  be  written 

1  =  t*2  Ne#[(y  +  7y)  erf  *  +J=  7  exp  <-y2)  *  7  |erf  ^(et)  - 

4“i7!exp(^2,!(lT)  •  ] 

If  we  now  consider  large  drift  velocities  then  we  need  the  case  for  large  y.  The 

2 

second  term  in  the  expression  rapidly  goes  to  zero  since  it  varies  as  exp  (-y  ).  The 

error  function  becomes  asymptotically  unity  for  large  values  of  the  argument,  thus 

the  constant  term  in  the  expression  becorr  r's  equal  to  y.  The  term  containing 

(eV/kT)  is  difficult  to  deal  with  since  it  contains  a.  However,  we  note  that  it  also 
2 

contains  exp  (-y  )  which  tends  to  zero  very  rapidly  for  large  y.  For  smaller  values 
of  (eV/kT),  where  the  second  and  higher  order  terms  can  certainly  be  neglected, 
we  have 


I  =  7tR2  Ne  Vg 


T 


1  eV 

7  ET 


Thus,  the  linearity  of  the  accelerating  characteristic  is  not  destroyed  by  large 
rocket  velocities.  In  fact,  the  linear  approximation  may  be  better  at  larger  rocket 
velocities,  since  the  linear  term  decreases  as  (erf  y)/y,  whereas,  the  square  term 
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2  2 

decreases  as  exp  (-y  ).  TTie  term  containing  V  also  has  a  in  its  coefficient,  and 
the  variation  of  a  with  y  has  not  been  calculated.  However,  it  seems  unlikely  that 
a  will  increase  with  sufficient  speed  with  y  to  offset  the  extremely  rapid  decrease 
due  to  the  exponential  factor.  Consequently  the  square  term  becomes  relatively 
less  important  in  comparison  with  the  linear  term. 

The  slope  (dl/dV)  decreases  as  (1  /v^)  and  thus,  in  the  limit,  becomes  zero. 
TTie  current  is  now  independent  of  the  voltage  and,  as  we  would  expect  from 
elementary  considerations,  has  the  value 

2 

IQ  =  irR  Ne  vQ 

Carrying  through  a  similar  expansion  for  the  retarding  case  gives 


I  =  ttR*  Ne 


erf  y  +  J 
V7T 


exp  ( -y2)  t  i  |erlr|(g)  - 


4  !exp  <-y2>!  (kt)  +  •] 

Note  that  this  expression  (valid  for  all  y)  tends  to  the  series  expansion  for 
exp  (eV/kT)  for  y  -*  0  as  it  should.  For  large  y  the  constant  term  and  the  term  in 
(eV/kT)  are  identical  to  those  for  the  accelerating  case  and  hence,  the  continuity 
of  1  and  (dl/dV)  at  V  =  0  (the  changeover  from  retarding  to  accelerating  )  is  pre¬ 
served  at  high  drift  velocities.  The  exponential  form,  however,  is  destroyed.  The 

2 

ratio  of  the  coefficients  of  V  and  V  in  the  expansion  is 


2 

1  e  y  exp  ( -y  ) 

ET  erf  y 

2 

As  v0  becomes  large  this  tends  to  zero  as  y  exp  (-y  )  and  the  second  order 
term  diminishes  very  rapidly  with  increasing  drift.  The  retarding  characteristic 
also  becomes  linear,  and  from  the  continuity  noted  above  the  straight  line  of  the 
accelerating  characteristic  is  extended  into  the  retarding  region.  For  a  fixed  value 
of  Vq,  if  the  voltage  is  increased  sufficiently,  x  will  eventually  dominate  over  y 
again  in  the  expression  and  the  curve  departs  from  linearity,  and  the  exponential 
form  is  regained. 
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10.  THE  MULTI  -GRID  SENSOR 


To  extend  the  analysis  for  a  case  where  the  incoming  particles  flow  from  the 
ambient  plasma  through  a  set  of  concentric  spherical  grids  to  a  central  collector, 
it  is  necessary  to  take  the  equation  for  I  (p.  11) 


I  =  2*/  Ne(f)  i  f  fv2  exp 


j-a(v  -  v0>2 


-  exp  | -a  (v  +  v 


,1] 


cos  c  sin  c  d£  dv. 


and  perform  the  final  integrations  according  to  .he  limits  on  v  and  £,  given  by 
energy  and  momentum  equations  in  transition  from  grid  to  grid. 

Consider  a  particle  at  grid  i  of  radius  r.  with  velocity  vi  at  an  angle  £.  to  the 
inward  radius  vector  arriving  at  grid  i  +  1  radius  r^  +  ^  (r^  +  1  <  r^  with  velocity 
v^  +  ^  at  an  angle  £.  +  y  If  potentials  on  the  two  grids  are  V.  and  V.  +  then  from 
energy  considerations 


i 


7  mv?  *  j  mvi+  2  +  e  (V.  +  ,  -  V.) 


1  2 

Defining  a  velocity  u^  by  the  relation  j  mui  =  eV^,  this  equation  becomes 


2  2  2  2 
visvi+it"i+rui  • 


From  angular  momentum  considerations, 


rivisui<r..ri+1v.  +  1sm<ri+1  . 


The  particle  will  just  enter  the  grid  i  +  1  if  it  arrives  at  the  grid  travelling 
tangentially,  that  is,  if  sin  +  1  <  1,  that  is,  if 


viri 

vi  +  i  avi—  smn  • 


which,  on  substituting  in  the  energy  relationship,  gives 


2  2  2  2  ri  .2 

vi  +  ui  '  ui  +  Is  vi - 2  SU1  «i 

ri+l 


24 


Hence,  for  a  given  value  of  v^  the  particle  passes  through  grid  i  +  1  if 


«in2  i+1 

sin  s 


ri 


1  + 


ui  "  ui  +  1 


vi 


We  note  that  the  above  relation  is  satisfied  for  all  £.  (o  <  if 


r. 


i  +  1 

"r? 

i 


u.  -  u. 


1  + 


i  +  1 


*  1 


or 


2,2  2  v 

vi<(u.  -ui  +  1) 


2 

ri+  1 


i  +  1 


9  6 

which  can  only  apply  for  uf  >  ui  +  1 ,  that  is,  accelerating  potentials  from  grid  i 

2  1  2  2  2 

to  grid  i  +  1  since  vi  s  0.  For  vi  >  (ui  -  ui  +  j)a .  +  j 


Limits  on  e.  are  then  given  by  the  limits  on  Vj  of  v^min)and  thus 


r 

- 

2 

2  2 

2 

2  2 

ri+ 1 

ui  -ui  +  , 

^  .2  ^  ri+i 

,  +  VV 

2 

1  2 

i-sm  ci~  ~J~ 

1  2 

ri 

vi(min) 

ri 

vi(max) 

«)  2 

For  retarding  potentials  ut"  <  u.  +  j,  a  lower  limit  is  imposed  on  vi  by  the 
condition  that  sin  c.  *  0,  hence 
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2  2  2 

iil  L-i  ,0 


2  >  2  2 
vi"ui+l-ui 


The  set  of  limits. 


u l  -  ui  +2ia  0 


(retarding) 


h  +  WH*. 


2  2  2 
0<  sin2c  <  1  +  — - 

1  vi 


u.  -ul+1<  o 


(accelerating) 


0<vffi  <u2-ui+X+I 


(ui-ui+2l»“i+l<  vfS< 


0  <  sin^e/S  1 


2  2  2 

0,sin2e.,r-Lti  l  +  M'jjLL 

‘  ri  [  vi 


on  v.  and  c imposed  by  the  condition  that  the  particle  reach  the  i  +  l11  grid,  can 
*  *  th 

be  extended  to  the  i  +  k  grid,  since  the  velocity  with  which  the  particle  emerges 
at  this  grid  is  a  function  only  of  v.,e.,  u.,  r.,  u.  .  ,  and  r.  .  . 

jt  11  1  1  1  »  K  1“t"K 

Hence  for  the  kin  grid 


2  _  2 

.2  <  1  i  +  k  -  ,  Ui  ui  >  k 

sm  ,  <  —r—  1  4  - j - 

ri  L  vi 


„2  r  2 

,  2  ,pl  .  ui 


tli  '2*  *  1  u  * 

For  i  =  0,  let  r^  =  r,  and  =  0.  Hence,  for  the  iin  grici  sin  f ^ — 2”  1  — 2 

dropping  the  O^1  subscript.  r  v 
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The  effects  of  the  various  grids  in  limiting  v  and  e  (the  initial  particle  parameters 
at  grid  0)  can  now  be  compared  and  the  integral  I  performed  according  to  the  way 
each  grid  (i  =  1,  2,  3,  etc. )  limits  the  values  of  v  and  e.  Since  the  lower  limit  on 

€.  is  always  zero,  the  e  integration  may  be  performed  prior  to  discussing  the 

1  2  2  2 
specific  applicable  limits  on  v,  writing  r  /r.  s^,  we  have 


Hence,  we  may  write  for  I, 


2  2 

The  limits  imposed  by  v  on  sin  e  can  be  illustrated  by  plotting  the  limiting 
curve 


9 

for  each  grid  of  radius  with  potential  Vi  (u.  =  2eVi/m). 

Various  potentials  will  be  considered  for  a  four-grid  case,  to  illustrate  the 
progression  of  limits  on  the  I  integral.  These  are  so  chosen  to  include  all  possible 
I  integrals  but  do  not  cover  all  possible  potential  differences  between  the  various 
grids.  The  resulting  equations  will  then  be  extended  to  cover  the  general  case  for 
N  grids. 

Consider  the  i^1  and  the  kth  grid,  where  k  >  i. 


Suw 


and 


.  2 
sin  e 
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The  intersection  point  of  these  two  curves  is  given  by 


v».v» 

k  l 


2 

R2  - 

P  i  -  ~y 
ri 


2  2 

.2  .2  ui  "  Uk 

sin  €  =  sin  e = - n - n - n - n 

V  ui  ■  pi  uk 


u2=!!Xi 

1  m 


and  for  the  intersection  point  to  lie  within  the  domain  of  the  integral  (e,  v)  then 


vika° 


and 


0<  sin  e.k<  1 


Since  ^  >  p\ >  1  tlien 


2  “k«»?  -» 

k 


Giving  the  restrictions  on  uf  of 


2  2 
uiauk 


for  0 


and 


for  real  v.^  and  where  i  <  k. 


2  2  Pi  '  1 

Ui  -  Uk  - 


K-1 


for  u^<  0 
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Before  turning  to  specific  cases  it  is  of  interest  to  note  that  for  the  outer  grid 
radius  r,  defined  to  be  at  zero  potential, 


a2  -  r  -  i 

P0  "T  1 
r 


Hence, 


and  u 


0 


=  0. 


From  which  it  can  be  seen  that,  for  the  outer  grid  to  contribute  to  the  I  integral, 
the  kl  grid  potential  must  be  such  at  u^<  0.  Hereafter,  the  grid  k  having  u^<  0, 

that  is,  eV.  <  0  or  the  sign  of  the  grid  potential  being  opposite  to  that  of  the  particle 

K  2 
charge,  will  be  termed  an  accelerating  grid.  Conversely,  one  having  u^  *  0,  will 

be  termed  a  retarding  grid. 

Since  the  above  equations  define  an  upper  limit  to  sin  e,  the  limits  of  integration 
lie  below  the  limiting  curves. 

Taking  a  configuration  of  four  grids  with  the  radii  having  the  ratios 
r2  :  rj  :  r|  :  r33  =■  1  :  0.  8  :  0.  4  :  0.  2, 


hence. 


2 

^\=~2  =  1.  25,  022  =  2.  5,  and  P*  =  5.  0. 

We  can  now  investigate  the  successive  limits  on  the  I  integral  by  applying  various 
combinations  of  voltages  to  grids  1,  2,  and  3. 

(a)  U3  >  u|  >  u2  >  0 

2  2  2 

The  limiting  curves  for  -  6,  Ug  =  4,  and  Uj  s  3  are  shown  in  Figure  5. 
Since  u2  >  u^,  no  curves  intersect  and  the  limits  on  v  go  from  point  A  to  infinity 
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(b)  uf>  u22>  u2  >  0 

2  2  2 

Figure  6  shows  curves  ior  uQ  =  1,  u0  =  2,  and  u.  =  3.  The  condition  for  real 

2  2  6  4  i 

intersection,  u.  ^  u,  ,  is  now  fulfilled  by  grids  1,  2,  and  3,  but  not  by  grid  0.  The 

1  K  2  2  2 
three  intersection  points  are  given  by  v^  =  4»  vj3=  H/3,  and  v2j  =  3  at  Bj,  B,  and 

B0  respectively. 

c  o 

The  limits  on  v9  now  commence  at  A  with  the  intersection  of  sin  e  -  0  and  the 
^  2 

most  retarding  grid  (in  this  case  grid  1)  intersecting  at  v  =3.  The  limiting  curve 

2  2  2 
is  then  given  by  Uj  =3  up  to  the  smaller  of  the  two  values  vJ2  and  v^2 .  Here, 

v13=  B  =  11/3  has  the  smallest  v^,  and  for  v^>  11/3,  the  limiting  curve  becomes 

n  2 

u£  =  1.  The  integral  I  is  now  the  sum  of  two  parts— from  A  to  B  along  Uj  =  3  and 
from  B  to«  along  u^  =  1,  and  grids  1  and  3  are  active  grids. 


E(v)  dv  + 


E(v)  dv 


where. 


rl  U1  -  r3  u3 


rl  ”  r3 


It  is  worth  noting  that  in  both  Eq.  1  and  Eq.  2,  where  the  active  grids  are 
retarding,  the  expression  for  I  is  independent  of  r,  the  outer  grid  radius. 
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2  2  2 

Figure  7.  Limits  of  Integration  for  0  >  >  u2  > 


(d)  0>  Ul2>  u22>  u2 

_  2  2  2  2  2 
Where  r  igure  8  now  shows  u3  =  -8,  u2  =  -6,  and  u1  =  -0.  25,  sin  <r  =  1  (uQ  =  0) 

is  the  limiting  curve  up  to  B,  the  smallest  intercept  of  vQ3  =  2(3^),  vQ2  =  4(B2> 

and  V021 =  1  (B),  making  grid  0  an  active  grid.  The  limit  is  then  given  by  u^  =  -0.  25 

up  to  C,  the  smallest  intercept  of  v12  =  5.  5  (C^  and  v^  =  7/3  (C).  The  final 

portion  is  from  C  to  ®  a'ong  u^  =  -8  making  grids  1  and  3  the  active  grids,  and 
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v 
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E(v)  dv  + 


2 

v 
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u3 

1-4 
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'01 


E(v)  dv 


2  v 
C1  T~ 


U1 


E(v)  dv  + 
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where 


u 


2 

1 


and 


2  2 
2  rl  U1 
f13 - 7 


u 


2 

3 


We  will  finally  take  two  examples  where  all  three  grids  become  active  and  the 
outer  grid  is  active  or  not  according  to  whether  all  active  grids  are  accelerating  or 
not. 
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2  2  2 

Figure  10.  Limits  of  Integration  for  0  >  Uj  > 


From  these  examples  it  can  be  seen  that  three  types  of  integrals  are  involved  in  the 
summation  for  I,  namely, 


The  summation  is  to  be  performed  over  the  active  grids  taken  in  sequence  from  the 
outermost  to  the  innermost  grids.  We  will  first  evaluate  these  three  integral  types 
and  then  develop  the  summation  for  I  for  a  general  case  having  n  grids. 

For  brevity  in  writing  out  the  resulting  equations,  let 

o 

exp  ( -ax  )  -  E(x) 
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and 


erf  (t£x)  =  F(x) 
where 


x  2 

erf(x)=2_  f  e_t  dt 

*  Jo 


by  definition. 


The  above  integrals  can  now  be  expressed  in  terms  of 


B  B 

IQ  =  J*  E(v  -  vQ)  d v  ~  j  E^v  +  vo^  ^ 


and 


B 

'  / 


I2  =  /  v2  E(v  -  vQ)  dv  -  f  v2  E(v  +  vQ)  dv 


/ 


which  are  both  standard  forms. 


E(v  +  vQ)  dv  *  ^  || j  [f(B  +  vQ)  -  F(A  +  vQ) 


B 

/ 


v  E(v  +  vQ)  dv  =  (B  -  vQ)  E(B  +  vq^  ‘  '  vo*  E^A  +  v0* 


-Mf)  (v0 +  4 )  [F(B  +  v  -  F<A  +  V 


Using  these  general  expressions  we  may,  upon  inserting  the  appropriate  limits 
A  and  B,  evaluate  the  three  desired  integrals  Ig,  I  ,  and  I  #  noting  that  for  1^, 
when  A  -►  0 
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E(A  +  v0)  -►  E(v0)  , 

and 

F<  A  +  vQ)  —  F(vq) 
For  I  when  B-*~®, 
E(B  +  vq)—  0  , 

and 

F(B  +  vQ)  —  1  . 


E(v)  dv  =  r2  r  T  *1 

*T 


B=v 
A  =  0 


Oi 


I  -  r2  j  V° 

h  T  2a 


E(vq)  -  E(voi  +  v0)  +  E(vq)  +  E(v0i  -  vQ) 


V01 

"2a"" 


E(v0i  -  vo'  -  E(v0i  +  V  ]  -Ui--(r)'  [F(v0i+v0> 

.'I! 


-  «*0i  -  vo>  -  F(vo)  -  F(v 


1 

=Tl(avn  1  2)(l)7  f<V-7  |:<V  - 


K  ♦*)/,? 


To—  r  F(voi  +V 


v0i+v0  ,  ,  vm  -  ■()  ,„  ,  ,  | 

—2a -  Mv0i  ’  v0)  4  —2a -  Mv0i  '  V  ' 


'  F<v0i  -  V 
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r2  v2 

r j  T 


f  Uf  1  ri2  f*k  2  rf  Uj  f  ]k 

|l  -  -ij-j  E(v)  dv  =  J  v  E(v)  dv  -  J  E(v)dv. 


2  t}  2  2 

^  pk 


F(vij +  V  -  F(vij  -  v0) 


+  1J2a  '  0  E(vij  -  v0> 


E(Vij  + 
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+  v0)  -  F(vjn 


Writing  out  the  resulting  expression  for  Example  5,  we  note  that  the  first 
integral  is  of  the  type  Im  with  =  Uj,  Vjk  =  v12;  the  second  has  v.j  =  v12,  v^k  =  v23; 

and  the  last  is  1^  with  Vjn  =  v23.  The  sequence  of  intercepts  is  Up  v12,  v23,  «  , 
giving  for  I: 


I 


12Za  °  E(v12  +  Vl  +-tT  °  Y*  [F(v23  +  v0>  '  F(v23  ’  Vl* 


.  v23 +  v0  ,,,  ,  v23  ’  v0  v 

+  — 2a h(v23  ’  v0> - 2a E(v23  +  V 


Collecting  terms  together. 


2((avo2  +  7-aui2)/M7r  . 

rlj  - 2a - '(«)  [F(ul  +  V0 


1  U1  +  v0 

)  -  F(Ul  -  vQ)  +  -  2a  E(Ul  -  vQ) 


^-2.  E(Uj  + v0)(  +  (r2  "  rl 


av2  +  i  -  av,2  \,  ,7 1 


F<v12  +  v0>  ' 


1  v12  +  vn  v12  -  v0 

-  F(v12  -  v0  >  +  -T, r —  E<v12  -  V  -  -TT1  E<v12  +  V  + 


+  (r 3  -  r2  )  j  (  — 27 - ~Vr)  [F(v23  +  v0)  *  J(v23  -  v0>  + 


'23  v0'  1  23  v0' 


.  v23+'o  c„  v  v23-v0  L 

+ - 2a E(v23  -  v0> - 2a  E(v23  +  V  • 


For  Case  6,  the  sequence  of  intercepts  is  0,  vQ1,  W2V  ®*  Sivin6»  in  the  same 

way  as  5 


,  P  [  2  f/av0  +  7\/jt\7  ,  ,  v0  ,  2  2*  f/av0  +  7  *  av01 

1  T  p  jl— 2a— X«  )  2F(v0>  H  27  2E(V0>  +  (rl  *  r  »  I - Za - 


Vrt1  +  vr 


[l)  F(v01  +  v0>  -  F(v01  -  v0)  J  +  012a  9  E(v01  •  v0)  -  012a  ° 


'01  ’0' 


E(v01  +  V  +  <r2  *  rl2)  |  • 


. .  |  +  (r  2  -  r  2) 


41 


p  =  £Ne  l  =  ?rNe  f8ET~  J_ 
Zajs  y  "  “2“  y  rnn  y 


42 


1 


where  P  =  2irNe 


Hence 


r5kf  1  ^  (r2 

7rm  y  /  ^  i 


'"■’if* 


+  y 


s(i-l),  i) 


■6 

7 


erfl*(t-i>.i+1'!-erfix(i-i).i*Y 


x(i-l),  f1 

- 7— 


expl(xu-l). 


exp  |  -(x 


(i-1). 
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Where,  if  the  first  active  grid  is  retarding  ,  that  is,  eV^—  0,  Uj  —  0,  summation 
starts  at  i  =  1  wi*h  r._j  =  Tq  =  0.  For  the  first  active  grid  accelerating,  eV^<  0, 
u2<  0,  summation  starts  at  i  =  0  with  r.  ^  =  r_j  =  0.  Note  that  with  'rQ  =  0 
(by  definition)  and  rQ  =  r  for  the  retarding  first  grid: 


and 


2  2  2  2 

2  r0  au0  -  rl  aul  2 

X01  2  2  aUl  ; 

ro  -ri 


and  for  an  accelerating  first  grid: 


,  2  2v  22 

r0  -  r-i  >  =  r0  =  r 


-10 


r-l  au-l  •  r0  au0 

- 2 - 2 - 

r-l  “  r0 


=  0 


and 
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2  2  2  2 
rQ  auQ  -  auj 


r0  ’rl 


The  method  used  to  determine  the  active  grid  sequence  i  =  0  or  1  through  n  from 
the  total  grid  sequence  K  =  0  through  N  may  be  summarized  as  follows: 

(a)  Compare  the  sequence  of  intercepts 


2  th 

k  =  1  through  N,  and  find  the  grid  K  =  kj  having  the  least  xqk*.  this  ki  grid  is  the 

2  th  L 

first  active  grid.  If  this  value  is  negative  or  zero,  this  kin  grid  is  retarding 
and  i  =  1.  If  two  or  more  grids  give  the  same  intercept,  the  one  with  the  higher  k 


is  chosen. 


(b)  Compare  the  sequence  of  intercepts 


2 

k  =  kj  +  1  through  N,  again  finding  the  grid  k  =  k2  having  the  least  this  grid 

is  now  the  second  active  grid  (i  =  2). 

(c)  Repeat  step  (b)  until  kn  =  N  when  the  active  grid  sequence  is  complete 
(i  =  n).  Note  that  the  N^1  (innermost)  grid  is  always  an  active  grid  because  its 

limiting  curve,  whether  accelerating  or  retarding,  tends  asymptotically  to 

2  2  2  2  2  2 
sin  c  =  1/0^  =  r^  /r  as  v  -►«>,  thus  having  the  least  value  of  sin  c  of  all  the 

grids  in  the  system  and  being  the  final  integral  in  the  sum  for  I. 

From  the  foregoing  discussion  of  a  system  of  concentric  spherical  grids  placed 

in  a  drifting  Maxwellian  plasma,  we  conclude  that  the  current  flow  through  the 

innermost  grid: 

(a)  is  a  function  of  the  plasma  parameters:  charged  particle  number 
density  N,  particle  temperature  T,  particle  mass  m,  particle  charge  c,  and  the 
drift  speed  vQ; 

(b)  is  not  necessarily  dependent  upon  all  grid  radii  and  the  potentials  thereon 
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(c)  if  one  or  more  grids  are  retarding,  is  independent  of  any  grids  of  radii 
greater  than  the  most  retarding  grid,  including  the  outermost.  If,  however,  no  grids 
are  retarding,  the  current  flow  is  dependent  on  the  outer  grid  radius  but  independent 
of  any  grids  of  radii  greater  than  that  of  the  "least  accelerating"  grid. 


II.  APPLICATIONS 

The  relations  derived  in  the  preceeding  sections  can  be  used  to  obtain  the  plasma 
parameters  from  the  observed  current -voltage  curve  for  the  probe. 

Since  the  various  expressions  all  involve  particles  of  only  one  sign  -  it  is  of 
interest,  at  this  point,  to  investigate  the  possibility  of  using  a  probe  which  responds 
to  particles  of  one  sign  only,  rejecting  the  other.  To  this  end,  it  seems  that  a 
reasonable  method  would  be  to  put  a  large  enough  bias  on  the  collecting  surface  to 
repel  all  particles  of  one  sign,  but  to  surround  this  element  with  a  perforated  grid 
so  that  the  effects  of  this  bias  will  not  penetrate  the  surrounding  plasma.  This,  then, 
gives  us  a  two -element  sensor.  The  final  conclusions  in  Section  10  stated  that  the 
inward -current  flow  is  a  summation  over  all  active  grids  and  always  includes  a 
term  containing  the  potential  and  radius  of  the  innermost  grid,  in  the  present  case, 
the  collector.  This  term  is  for  an  attractive  potential  V., 

T2  =  (rl2  •  r22)|(?+  y2  '  xl2)'^[erf!x12  *  A  -  erflx12  -  dj  + 

+  ~z — ~~exp|-(x12  -  y)2|-  122  exp  j-(x12  +  r>2  jJ 

2  2  2 

It  can  be  shown  that  if  x^2  »  V  and  xi2  >:>  *  the  above  term  tends  to  zero, 
and  the  current  flow  will  be  independent  of  the  inner  grid  potential  for  the  attracted 
sign  particle.  For  the  retarded  particle,  however,  the  total  current  flow  has  only 
one  term  (under  the  assumption  that  |v2(  »  |Vj|},  namely. 


J 8kt*  1_ 
1  irm  y 


-  erf 


x02  +  V  ( 

2 -  exp  r(x02  •  Y) 


<x02  + 


+ 


2  2  2 

Again,  if  xQ2  »  y  and  xy2  »  1,  then  the  total  current  will  tend  to  zero  and  there 
will  be  no  contribution  from  the  retarded  sign  particle. 

Putting 
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On  the  assumption  that  the  bias  potential  on  the  collecting  electrode  satisfies 
the  set  of  conditions,  we  may  proceed  with  the  analysis  on  the  basis  of  having  a 
spherical  probe  of  radius  R(=  rj)  at  a  potential  V  (=Vj)  measuring  particles  of  one 
sign  only. 

Not  all  of  the  particles  which  are  incident  upon  the  sphere  with  radius  equal  to  that 
of  the  outer  (perforated)  grid  will  contribute  to  the  current  since  some  will  strike 
the  grid  wires.  The  observed  current  is  reduced  by  the  transmission  factor  of  the 
outer  grid  t,  .  The  true  current  I  will,  therefore,  be  the  measured  current  divided 
by  4  .  A  further  complication  arises  from  the  fact  that  all  potentials  in  the  preceding 
sections  are  referred  to  the  plasma  potential  which  is  unknown.  In  practice,  the 
sweep  potential  applied  to  the  probe  must  be  referred  to  some  other  point  in  the 
system  such  as  one  of  the  electrodes  in  a  discharge  tube  or  the  skin  of  a  rocket  in 
a  rocket  experiment.  The  plasma  potential  is  one  of  the  unknowns  we  would  like 
to  determine.  In  the  following,  V  is  understood  to  be  measured  from  the  reference 
potential,  and  it  is  assumed  that  the  potential  difference  between  the  reference 
electrode  and  the  plasma  remains  constant  regardless  of  the  potential  applied  to  the 
probe. 

First,  we  consider  the  case  for  no -drift.  In  the  retarding  region,  for  the 
relation  between  I  and  V,  we  can  write 

1  =  IBexp  (’It)  • 

To  simplify  the  discussion  we  will  consider  the  case  of  electrons,  and  defining 
a  =  |  e  |  /kT,  we  have 

I  =  IB  exp  (a  V) 

In  the  "near  accelerating"  region,  where  the  current  may  be  represented  by 
the  first  two  terms  of  the  Taylor  expansion,  we  can  write 

I  =  a(V  -  VA)  . 

We  have  now  introduced  four  constants:  Ig,  a,  a,  and  V^.  They  are  not, 
however,  all  independent,  since  we  have  shown  that  the  gradient  (dl/dV)  must  be 
the  same  at  the  plasma  potential  for  both  retarding  and  accelerating  expressions. 
Differentiating  the  above  expressions  gives 


retarding 
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^  =  a  accelerating 

At  the  plasma  potential  Vq  the  current  flow  will  be  Iq.  Equating  these  two 
gradients  at  the  plasma  potential  gives 

aI0  =  a' 
or 


This  gives  us  a  means  of  obtaining  the  current  at  plasma  potential  Iq.  We  plot 
log  I  versus  V,  which  will  be  a  straight  line  in  the  retarding  region.  The  slope  of 
this  line  is  a.  A  plot  of  I  versus  V  will  give  a  straight  line  in  the  near  accelerating 
region.  If  the  conditions  for  the  particular  experiment  are  appropriate,  there  will 
be  a  reasonably  large  range  of  voltage  over  which  the  two-term  Taylor  expression 
will  be  a  good  approximation,  and  this  straight  line  will  be  readily  located.  Its 
slope  is  a  and  the  current  at  plasma  potential  is  obtained  as  the  ratio  of  the  two 
slopes.  Note  that  the  determination  of  this  current  did  not  involve  an  accurate 
determination  of  the  plasma  potential.  The  slope  a  is  simply  related  to  the 
temperature,  hence,  electron  temperature  is  obtained  from  the  equation 


The  electron  density  is  given  by 
IQ  =  2.  48821  X  10‘16  (4?rR2)  Nn/t 


The  intercept  of  the  straight  line  can  also  be  determined,  and  from  the 
equation 

I0  =  a(V0-VA)  , 
we  have 


which  gives  the  plasma  potential. 
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We  have  now  determined  three  plasma  parameters  from  the  current-voltage 
characteristic  of  the  probe,  namely,  the  electron  temperature,  the  electron 
density,  and  the  plasma  potential.  'Riis  is  the  maximum  number  that  can  be  obtained. 
TTie  curve  fitting  process  introduces  four  parameters  as  already  noted.  However, 
there  is  one  relationship  between  them,  that  is,  the  condition  to  match  the  slopes  of 
the  two  curves  at  plasma  potential.  Thus,  there  ar*»  only  three  independent  parameters. 

This  particular  point  has  more  significance  for  positive  ions  when  the  mass  m 
is  also  an  unknown.  In  this  case  only  three  of  the  four  unknowns  can  be  obtained 
from  the  probe  characteristic;  the  temperature  is  obtained  from  or,  the  plasma 
potential  from  a  and  V^,  and  then  using  a  in  conjunction  with  a  to  give  Iq,  we  get 
N/N/m.  Hence,  either  m  must  be  known  so  that  N  can  be  calculated,  or  vice  versa. 

When  the  drift  velocity  cannot  be  neglected,  the  equations  are  much  more 
involved.  The  straight  line  part  of  the  accelerating  characteristic  now  extends  into 
the  retarding  region  as  noted  earlier  but  even  with  large  drift  velocities,  if  x 
(that  is,  voltage)  is  made  large  enough,  a  curvature  in  the  characteristic  will 
eventually  occur.  This  is  likely  to  be  of  limited  practical  importance,  since 
making  the  voltage  large  for  the  retarding  case  implies  very  small  currents  which 
will  probably  be  less  than  the  noise  level.  Even  if  this  can  be  observed,  there  is 
still  no  simple  analytic  procedure  which  can  be  apolied  to  it  comparable  to  that 
which  we  can  use  in  the  case  for  no -drift.  However,  there  is  still  the  straight  line 
part  of  the  I-V  characteristic  (which  now  extends  over  parts  of  both  the  retarding 
and  accelerating  regions)  from  which  a  slope  and  an  intercept  may  be  obtained. 

The  slope  obtained  from  this  line  is  given  by 

2xR2  Ne2  f  /  JnT  \ 
mv0  ^  0  1ZFT  ). 

This  expression  contains  the  density  N,  the  mass  m,  and  the  temperature  T, 
all  of  which  are  unknowns.  Note  that  the  plasma  potential  does  not  enter  into  this 
expression.  To  obtain  an  expression  for  the  intercept  suppose  we  fit  the  expression 

I  =  a(V  -  VA) 

to  the  data  so  that  is  the  intercept  on  the  voltage  axis.  For  the  plasma 
potential  we  have 

‘o  5  a<vo  -  V  • 


and  since  a  is  the  slope, 
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Using  the  expressions  already  developed  for  IQ  and  (dI/dV)Q  we  obtain 


where  the  upper  sign  applies  to  positive  particles  and  the  lower  to  negative  particles. 
This  expression  involves  the  three  unknowns:  plasma  potential  VQ,  the  mass  m, 
and  the  temperature  T,  which  occur  in  y.  Note  that  the  density  N  does  not  enter 
into  this  expression. 

We  now  have  two  relations  between  four  unknowns.  These  enter  into  the 
expressions  in  such  a  way,  that  if  independent  estimates  can  be  made  of  any  two 
unknowns,  the  other  two  can  be  determined.  However,  some  cases  will  be  much 
simpler  to  compute  than  others.  For  example,  solving  for  VQ  and  N,  when  m  and 
T  are  known,  is  straightforward.  On  the  other  hand,  solving  for  m  requires  a 
numerical  iteration  process. 

As  an  example  of  how  the  procedure  might  apply  in  practice,  consider  the  case 
of  a  rocket  carrying  two  probes  — one  measuring  electrons,  and  the  other,  positive 
ions.  In  practice,  due  to  the  high  thermal  velocity  of  the  electrons,  the  rocket 
velocity  can  be  neglected  for  the  electron  probe,  however,  for  the  ion  probe,  y  is 
likely  to  have  a  value  of  about  unity.  For  the  electron  probe  there  will  be  only 
three  unknowns  since  the  electron  mass  is  known.  From  the  retarding  part  of  the 
characteristic  we  obtain  electron  temperature,  and  using  the  accelerating 
characteristic  as  well,  gives  density  and  plasma  potential.  We  might  now  assume 
electrical  neutrality  so  that  the  ion  density  is  equal  to  the  electron  density.  Since 
the  plasma  potential  has  already  been  determined,  we  have  VQ  and  N  for  the  positive 
probe  and  can  use  the  probe  results  to  obtain  mass  and  temperature.  Alternatively, 
we  might  assume  that  the  ion  temperature  is  equal  to  the  neutral  temperature  which 
is  known  from  independent  experiments  and  determine  N  and  m  from  the  probe 
data,  or  we  might  use  independent  results  on  composition  to  obtain  m  and  determine 
N  and  T. 


— rmbri*VV~&\ 


r*  =*->  -r*B*gg 
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Appendix 


We  list  below  the  equations  we  have  derived,  together  with  some  of  their  special 
forms.  The  following  notations  are  used  throughout: 

I  s  current  flowing  to  the  spherical  collector, 

V  =  potential  difference  between  collector  and  plasma, 

R  =  radius  of  collector, 
r  =  radius  of  sheath, 

N  =  charged  particle  density  in  the  plasma, 
e  =  charge  on  the  particle, 
m  =  mass  of  charged  particle, 
k  =  Boltzmann  h  constant, 

T  s  Temperature  of  charged  particle  being  collected, 


v0  *  drift  velocity  of  plasma  relative  to  collector, 
a  =  m/2kT, 

r  =^v0 

*=  R2/(r2  -  R2), 

I|eV/kT  for  retarding  voltages, 

(R2 /(r2  -  R2)(|eV|/kT)  for  accelerating  voltages, 


erf  y  = 


2 


/ 


exp(-t^)dt. 


i 

i 

! 


! 

I 

; 


i 

i 

;■ 

l 
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(a)  Retarding  Express  ion:  - 


I  =  irR  Ne 


•jnr  7  [?  (?  +  ?2  ‘ x2) ! erf  <x  +  Y) ' erf  (x  ‘  Y,l+ 

+  J  (x  +  y)  exp  |-(x  -  y)2J  -  j  (x  -  y)  exp  J-(x  +  y)2 


(b)  Accelerating  Expression:- 


^  (1  +  2y2)  erf  y  +  yexp  <-y2)  -  j  £(  £  +  y2  -  x2) 

erf  (x  +  y)  -  erf  (x  -  y)j  +  ^  (x  +  y)  exp  -(x  -  y)2j- 


1  2 
-  j  (x  -  y)  exp  -(x  +  y) 


II] 


(c)  Taylor  Expansion  of  Accelerating  Expression: - 


I  =  ttR  Ne  v 


0  \(l  +  ^r)erf  (Y)  +Jr  h  exp  ('y2)  *  7  I erf  <Y)|  (kt)  ' 


’?  -I  ^  |exp  (~y2>{ (fr)  +  •••  |  . 


where  the  upper  sign  is  for  positive  particles,  and  the  lower  for  negative. 


(d)  Current  at  Plasma  Potential  (V  =  0):- 


I0  =  irR2  Ne  vQ  [(l  +  -iy)  erf  y  +  J.  i  exp  (-y2) 


(e)  Derivative  of  tho  Current  at  Plasma  Potential: - 


/dl  \  - 

WJ  „  " 


2ttR2  Ne2 
mvn 


erf  y  . 
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where  the  upper  sign  is  for  positive  particles,  and  the  lower  for  negative. 
When  the  drift  velocity  is  zero,  the  above  expressions  reduce  as  follows: 

(f)  Retarding  Expression:- 


=  4ri?2Ne  VHji  exp  ('  ct) 


(g)  Accelerating  Express  ion:  - 


I  =  4jrr  Ne 


VJT 


2  -  R2 

1 - ^  exp 


/  R2 

l  7~Tr 


eV 

T  FT 


(h)  Taylor  Expansion  of  Accelerating  Expression:  - 


I  =  47rR  Ne 


-fW  [lTeV  1  R2  i 

[1T  FT  T  .  r2  \] 


eVV 

FT) 


(i)  Current  at  Plasma  Potential:- 


V4ril2Ne 


Note  that  the  average  particle  velocity  in  a  Maxwell  distribution  is  given  by 


so  that  the  above  current  can  be  expressed  as 
IQ  =  7rR2  Ne  v  . 


This  checks  with  the  result  obtained  from  elementary  kinetic  theory—  the  number 
of  particles  crossing  unit  area  in  unit  time  is  given  by  Nv/4. 

The  expression  for  In  contains  the  universal  constants  e  and  k  and  the  numerical 

u  _.U  -23 

values  may  be  inserted.  Taking  e  =  1.  60210  X  10  coulomb,  k  =  1.  38054  X  10 

-31  * 

joule/°K  and  considering  electrons  for  which  m  =  9.  1091  X  10  kg,  we  obtain 
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IQ  =  2.48821  X  10"16  (4ttR2)  N 

where  Iq  is  the  current  in  amperes,  R  is  the  distance  in  meters,  and  N  equals  the 
particles  per  cubic  meter. 

Considering  positive  ions  of  atomic  mass  M,  and  taking  the  atomic  mass  unit  to 
be  1.  65979  X  10~27kg,  we  have 

IQ  =  5.  82906  X  10"18  (4;rR2)  N  ^  . 

(j)  Derivative  of  Current  at  Plasms  Potential:- 


(w)  =  T  4,rR2  N®2  ' 

0 

Substituting  numerical  values  as  above,  we  have  for  electrons 


=  2.  88753  X  10"12  (4ttR2)  £1 


and  for  positive  ions, 


=  6.  76453  X  10"14  (4ttR2)  — 
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